Non-thermal images of SN 1006: from radio to 7-rays 



Vasyl Beshley 1 

Abstract: Supernova remnants (SNRs) are believed to be the main sources of galactic cosmic rays. Discovery of the non-thermal 
component in X-ray spectrum of SN 1006 in 1995 and detection of a number of SNRs by H.E.S.S. strengthen the investigation of 
SNRs. SN 1006 remains to be one of the most interesting objects for high-energy astrophysics. Electrons accelerated by the shock 
are the source of the non-thermal radiation in radio, X-rays (via synchrotron emission) and 7-rays (via inverse-Compton process). 
Experimental images of SN 1006 are known in all these bands, including the very-high energy 7-ray range. An important task is 
therefore to model the distribution of the surface brightness in SNRs. We develop a method for synthesis of SNR maps due to the 
non-thermal radiation of electrons in radio, X-rays and 7-rays. In particular, the method takes into account the injection of particles 
and the behaviour of magnetic field at shocks with different obliquities as well as the radiation losses of electrons downstream of the 
shock. The method is used to model images of SN 1006. The images well correlate with the observations in radio and X-ray but not 
in 7-ray range. For coincide in 7-ray range it is needed to increases parameters rj and Br. 
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1. Introduction 

It is common believe that SNRs are the main sources 
of cosmic rays at least up to the knee at 3 x 10 15 eV and 
even beyond jBerezhko & Volk||2006) , jBlandford & Eich- 
|ler| [l987|. Electrons accelerated at the shock of SNRs ra- 
diate in wide range from radio to 7-rays. The SN1006 is 
the prime example of such an object. This object is a type 
of la SNR located at a distance of 2.2 kpc (Winkler et al. 
[2003) . 

The remnant of SN 1006 was first identified in the 
radio jGurdner & Miln| |1965) . In 1982, the first paper 
went describing were published about describe the first 
radio images of SN 1006 and Tycho (Reynolds & Chava-| 
|Iier)|l982) . 

Synchrotron emission from accelerated electrons is de- 
tected up to the X-rays in SN 1006 jRothenflug et al.|[2004| . 
In 1995, the first papers went out about domination of 
non-thermal X-ray component in some regions of SN 1006 
jKoyama et"aT][li95) . 

Investigation of the high energy 7-rays was the next 
step. The first observations in this range are done by 
the CANGAROO telescope jTanimori et~aT) [1998). The 
first images of SN 1006 as a TeV source was reported by 
H.E.S.S. team in 2008. 

It is suitable to model this object because it has the 
simple symmetrical distribution of surface brightness in 
all ranges. Almost spherical form of SN 1006 provide ev- 
idences that it develops in a uniform interstellar medium 
(ISM). Therefore, we model SN 1006 and synthesize it is 
images in different bands as SNR in uniform ISM and uni- 
form interstellar magnetic field (ISMF). 



2. Description of The Model 

SN 1006 is on the adiabatic stage of evolution. There- 
fore, for description of hydrodynamics we use the Sedov 
self-similar solutions (Sedov 19591. The calculations of the 
magnetic field in the SNR interier and on the shock, are 
bathed on the model of Reynolds 1 19981. We use the classic 



model of kinetics of the charged particles. The spectrum 
of relativistic electrons is power-law with the exponential 
cut-off: 



N(E)dE = KE- S exp (-E/E„ 



)dE 



(1) 



where s is a spectral index and E max the maximum en- 
ergy of electrons. Distribution of spectrum of electrons is 
downstream of the shock is described in pa per |Reynolds| 
1998) and generalized in jPetruk & Beshley[[2008) ~ 

It is also needed to take into account the azimuth dis- 
tribution of the maximal energy. It is not isotropic as it 
is given by observations jRothenflug et aL 2004) . The 
Reynolds developed three models of maximum energy: 
l)by the electron radiative losses (due to synchrotron and 
inverse Compton, loss-limited), 2) by the limited time of 
the acceleration (time-limited) and 3) by properties of mi- 
crophysics (escape-limited). For description of distribu- 
tion of maximum energy we use the second model (time- 
limited) with parameter 17, defined as A = f]r g , where A is 
the mean free path and r^-gyroradius and r\ the "gyrofac- 
tor". 

[Reynolds] j!998) considered thee models for in- 
jection: quasi-parallel (qpar), isotropic (iso), quasi- 
perpendicular (qper). Morphology of SNR under quasi- 
perpendicular injection is similar to morphology for 
isotropic case {Orlando 2007). Ex perimental data reject 
quasi-perpendicular model jPetruk 2008), if the SNR is as- 
sumed to be in uniform ISM. Therefore we consider only 
isotropic injection. 

Spectral index. With the use of radio observation it is 
possible to find the value of spectral index. Spectral index 
in this case is s 



2.2 jRothenflug et al. 
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For the account of radiation losses of 
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Radiative losses 
electrons we use the parameter Ef, Ef = 637 (B 1 tE 
is the reduced fiducial energy, where t is the age of SNR, 
B the magnetic field strength after parallel shock. This pa- 
rameter is a measure of importable of radiative losses in 
modification of the electron spectrum (Reynolds, 19981. If 
< Ee < 1 then the radiative losses are very important if 
Ef > 1 radiative losses may be neglected. In our calcula- 
tions Ef = 1. 

Other parameter is orientations of ISMF in respect to 
the observer {(p ). Aspect angel equals 70° jPetruk 2008) . 
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In the present paper we consider two mechanisms of 4. Conclusions 
radiation of electrons. This is synchrotron (Petruk & Besh 
ley 1 |2007), {Petruk & Beshley| |2008) end inverse Comp- 



ton(IC) {Petruk 
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The surface brightness at an arbitrary point of SNR 
is calculated as integral along the line of sight, taking into 
account distributions of the magnetic field, density of elec- 
trons etc. 

3. Results 

Using results of observation we impose limitation on 
the parameters of model. 

Maximum energy. Distribution of the maximum en- 
ergy is taken from the observations of satellite XMM- 
Newton. For different azimuthal angle, E max is differ- 
ent. Observations agree well with time-limited model of 
|Reynolds|{i99"8) with rj = 1. 

Using the classic model of kinetics of electrons the 
maps of surface brightness distribution were built for 
three ranges: radio, X-ray and 7-ray. The model and ob- 
servational images of SN 1006 are shown on Fig. [I] The 
limbs of brightness in the simulated images coincide in 
radio and X-ray ranges. However we have no coincidence 
in 7-rays. The possible of such behaviour are discussed 
in JPetruk et al.J |2009) . The location of limbs, besides of 
azimuthal variations of injection efficiently and magnetic 
field compression, depends an a number of factors: 

• azimuth distribution of the maximum energy of elec- 
trons given by parameter r/ 

• radiative losses controlled by the parameter Ef 

• non-uniformity of the interstellar medium 

• non-uniformity of the magnetic field. 

With the growth of parameters t] and Ef, the maximal 
energy increases with the azimuthal angle and the radi- 
ation losses decrease. We expect that the values of these 
parameters t] and Ef must change in a range 5-10 and 30- 
60 respectively, in order to provide coincidence of limbs 
also in 7-rays. 



We considered the two models of emission of elec- 
trons: synchrotron (in radio, X-rays) and IC (in 7-rays) for 
the calculation of maps of surface brightness. Images are 
built in these three ranges (radio, X-ray and 7-ray) for SN 
1006. The observational and theoretical images in X-rays 
and radio coincide well, but in 7-rays does not correlate. 
To match the images we should account for the following 
modifications in our model: 

• introduce non-uniform of ISM and ISMF 

• for consider rj in the range 30 < rj < 60 

• do such the role of smaller radiative losses {Ef > 30). 
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